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Abstract 

A characterization of positive operators in finite dimensional complex vector spaces based on the 
Routh-Hurwitz criterion. 
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1 Introduction 

This article gives a characterization of the positive operators within the class of self-adjoint operators in 
finite dimensional complex vector spaces. The result is based on the Routh-Hurwitz theorem [ 1 1 which gives 
a criterion to determine if all (complex) roots of a real polynomial have negative real partQ 

Besides to provide a test for the positiveness of self-adjoint operators (in finite dimensional complex 
vector spaces), Theorem Q~2] also makes possible to deal with problems where positive operators occur as 
variables or unknown parameters. Exactly this situation was the motivation for this work. The problem 
had raised in a research about quantum dynamics, where the aim was to find conditions which guarantee 
the conservation of energy in dissipative quantum systems [2|. Since I could not find any reference to the 
desired result, I'm providing one! 

The article' structure is simple: section previews basic facts about positive operators and characteristic 
polynomials; section [3] presents the Routh-Hurwitz criterion and an analogous result which determines if 
all (complex) roots of a real polynomial have non-negative real part. In the section [4] the main theorem is 

'The Routh-Hurwitz criterion is well known and used to characterize stable dynamical systems. 



1 



stated and demonstrated, combining facts presented previously. Section [5] develops the application of the 
main theorem in two and three dimensions. 



2 Positive operators 

Let (J^, (, }) be a Hilbert space, where (, } is a inner product (anti-linear in the first entry and linear in the 
second entry). || || denotes the norm associated to this inner product: 



||9ll = VW> , V<p€JT 

Here, an 'operator in M" means a bounded linear operator from Jrf? to M '. (After some generalities, Jtf 
will be specialized to C", with n being a positive integer number.) 
In definitions below, I follow the references [3 |, [4| and |5|. 

Definition 1 (adjoint, self-adjoint operator, positive operator) Let A be an operator in Jtf. 
The adjoint operator of A is the operator A* : — > Jff such that^ 

(A>,y/) = (<p,Av/) , V<p,v/e 

A is self-adjoint when A* — A, i.e., 

(<p,Ay/) = (A(p,y/) , V<p,y/-£JT 

A is positive when 

(<p,a<p) >o , v<p Gje 

Proposition 2 Any positive operator is self-adjoint also. 

This proposition follows from polarization identity [3 p. 171 -172] J3] p.4], which gives the inner product 
in terms of the norm|l 

(V,<P) = 1{(\\V+<P\\ 2 -\\V-<P\\ 2 ) -'(llr + '<Pl| 2 -||V'-'<Pl| 2 )} , Vv/,<pe^ 

My aim is to characterize the positive operators within the set of self-adjoint operators. The argument of 
Section@]uses some definitions and results of spectral theory, of which I recall what is relevant. 



Definition 3 (eigenvalue, eigenvector, characteristic polinomial) Let A be an operator in An 'eigen- 
value of A ' is a number A G C such that ker [Xl — A) =^ {0}; in this case, the elements ofker (Xl — A) \ {0} 
are called 'eigenvectors of A associated to the eigenvalue A'. Explicitly: <p E Jt? is an eigenvector of A 
associated to the eigenvalue A E C if, and only if, 

<p / and Acp = A<p 

The 'spectrum of A' is the set of its eigenvalues: 

Oa '■= {eigenvalues of A} 
2 The existence of the operator A* follows from Riesz. representation theorem [ 4 p. 39] [5] p. 12, p. 3 1]. 

3 I note that positivity would not imply self-adjointness if the underlying field were the real numbers, because inner producs in real 
vector spaces cannot be writen in terms of the corresponding norms. 
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When J4? has finite dimension, the 'characteristic polynomial of A ' is defined by 

p A (z):=dct(zI-A) (1) 

In this situation, the eigenvalues of A are exactly the roots of the characteristic polynomial of A: 

VA G C.\ker(A/-A) ^ {0} <=^ p A {X)=0 

Proposition 4 Assume that Jf? has finite dimension and let A be an operator in Jff. 

If A is self-adjoint, then its characteristic polynomial is real (i.e., its coefficients are real numbers). 

Proof. Let p A (z) be the polynomial whose coefficients are the corresponding complex conjugated of the 

coefficients of the characteristic polynomial p A (z)," then: 



Pa (z) = Pa (z) , Vz G C 

Using the algebraic properties of the conjugation * (which assigns each operator to its adjoin), it follows: 



p A {z) = 6et{zI-A)=det(zI-A)* = det(zI-A)=p A (z) , Vz G C 
Combining the two identities above, it follows: 

PA (x) = p A (x) , Vx G K 

Therefore, the polynomials p A (z) and p A (z) have the same restriction to R; this implies that they have the 
same coefficients, from what follows that the coefficients of p A (z) are real numbers. ■ 

Proposition 5 Let A be an operator in the Hilbert space Jfif. 

i) If A is self-adjoint, then its eigenvalues are real numbers; 

ii) If A is positive, then its eigenvalues are non-negative real numbers. 

Proof. Let A be an eigenvalue of A and let <p G Jti? ((p ^0) an eigenvector of A associated to A. 
If A is self-adjoint, then 

x _ {<PM _ {(PA<P) _ (A(p,(p) _ {k(p,<p) _- x 
\\<P\\ 2 \\<P\\ 2 \\<P\\ 2 \\<P\\ 2 



If A is positive, then 



x = (jMj>) = (<P,A<P) >Q 
II«P|| 2 ll«Pl| 2 " 



Theorem 6 (Espectral) Assume that J%? is separable and let A be a (bounded) operator in Jrff. 
Then, A is self-adjoint if, and only if, there exists a sequence of orthogonal projections (P\,P2,...) 

P k Pj = PkPj = 8kjPk ,W^j = 1,2, ... 
and a sequence of real numbers (we can suppose to be distinct) (Ai,A2, ...) such tha^ 

1=^ , A=^hp k 

k k 



When the dimension of M J is infinite, the sum can be a serie which converges in the norm topology. 
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In this case 



a A = {A 1; A 2 ,...} 

and 

Im(P jfc )=ker(A jfe 7-A) ,VJt=l,2,... 

Theorem 7 Assume ^ is separable and let Abe a self-adjoint operator in J^f. 
Then, A is positive if, and only if, its eigenvalues are non-negatives. 

Proof. The implication was already demonstrated in the Proposition^ so, here I prove the reciprocal sen- 
tence (even though it also follows from the identity below). Take the spectral decomposition of A according 
with Theorem® then, for all (p G Jif it holds 




Therefore, if all eigenvalues of A are non-negatives, the last sum above is non-negative also. This proves 
that A is positive, according with Definition^ ■ 

From this theorem, we see that one can obtain a characterization of the positive operators within the class 
of self-adjoint operators in C" using any characterization of the real polynomials whose complex roots have 
non-negative real part. This points to the Routh-Hurwitz criterion. 



3 Routh-Hurwitz criterion 

Let a polynomial with complex coefficients, 

P (z)=b Q z" + b l z"- 1 
The Hurwitz determinants of p (z) are defined by 



-b n -iz + b„ 



At := det 
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where 



bj = , V; eN, j>n 



In particular, for n — 3 this definition means 



k = 1,2, ...,n 



b b \ , bi b 3 

Ai = b\ , Ao = det I , 1 , 3 ) , A 3 = det b b 2 
b ° b * ) \ b\ b 3 



(2) 



(3) 
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Theorem 8 (Routh-Hurwitz criterion) Consider the polynomial (0 has real coefficients and that its lead- 
ing coefficient is positive (bo > 0). Then, the roots of this polynomial have negative real part if, and only if, 
the Hurwitz determinants of the polynomial are positive: 

A k >0 , Vfc= l,...,n 

For a proof of this theorem I refer to |6, p. 231]. 

In the sequel, I will state and demonstrate theorems analogous to the Routh-Hurwitz criterion of Theorem 

m 

By convention, the phrase "w G C is a root of the polynomial p (z) with multiplicity zero" means that "w 
is not a root of p {£)"■ 

Proposition 9 (extended Routh-Hurwitz criterion) Consider the polynomial (O has real coefficients and 
that its leading coefficient is positive (bo > 0). Assume that the number zero is a root of p (z) with multi- 
plicity no, < «o < n. Then, the non-zero roots of p (z) have negative real part if, and only if, the Hurwitz 
determinants of the polynomial with indexes from 1 to n — no are positive: 

At > , Vk = 1 ,...,« — «o 

Proof. The case no — reduces to the Routh-Hurwitz criterion, Theorem |S} therefore, I will assume that 
no > 1. 

From the hypothesis, the coefficients of the monomials of p (z) having degree less than no are equal to zero 

b„-„ 0+l =...=b n =0 (4) 

/ define the auxiliary polynomial 

p(z):= boz n - na +b lZ n - n °- 1 + ... + b n - no - lZ + b n - no 

Therefore 

p(z)=z no p(z) 

The polynomial p (z) has degree n — no and its roots are equal to the non-zero roots of p (z); further, the 
Hurwitz determinants of p (z) and p (z) are also equal due to @ and Definition [3] As happens with p (z), 
p (z) is a real polynomial with positive leading coefficient; therefore, the thesis follows from Theorem [S] 
applied to p (z). ■ 

Remark 10 In the case p (z) = z n ( the number zero is a root of the polynomial with multiplicity equal to the 
polynomial's degree), all Hurwitz coefficients of p (z) are zero and (evidently) the polynomial has not any 
non-zero root. 

Proposition 11 (symmetric of the extended Routh-Hurwitz criterion) Consider the polynomial (0 has 
real coefficients and that its leading coefficient is positive (bo > 0). Assume that the number zero is a root of 
p (z) with multiplicity no, < «o < «• Then, the non-zero roots of p (z) have positive real part if, and only if, 
the Hurwitz determinants satisfy the condition: 

(_l) 1 +LV2J At>0 , vJt=l,2,...,n-no (5) 

where \k/2\ is the greatest integer less then or equal to k/2: 

[x\ := max{m 6 Z; m <x) , Vx 6 K 
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Consequently, the conditions @ are necessary and sufficient to guarantee that all roots of the polynomial 
p (z) have non-negative real part. 
Proof. I define the auxiliary polynomial 

p(z):=(-lTp(-z) 

Then, the non-zero roots of p (z) have positive real part if, and only if, all non-zero roots of p (z) have 
negative real part. The coefficients of p (z) are given in terms of the coefficients of p (z) by 

bk = (-1) , Vfc = 1,2,..., 7i 

Therefore, the Hurwitz determinants of p (z) and p (z) satisfy the identified 

A fc =(-l) 1+L * /2J A* , V*=l,2,...,fi 

As happens with p (z), p (z) is a real polynomial with positive leading coefficient (bo — bo) and the number 
zero is a root of it with multiplicity no; therefore, the thesis follows from the Theorem\9\applied to p (z). ■ 



4 Characterization of positive operators 

Theorem 12 (characterization of positive operators) Let Abe a self-adjoint operator in C and let no £ 
{0, 1, ...,n} be the multiplicity of the number zero as a root of the characteristic polynomial of A. Then, A is 
a positive operator if, and only if, 

(_l)i+L*/2J Ajfe>0 , V£=l,2,...,n-rt 

Proof. The result follows directly from the combination of the Theorem\7\and Theorem\ll\ ■ 

Naturally, this result motivates us to know some formula for the coefficients of the characteristic poly- 
nomial of operators. So, I highlight the following results: 

Theorem 13 (coeflcients of the characteristic polinomial) Let Abe a self-adjoint operator in C" and let 
its characteristic polynomial be 

p A (z) = det (zl - A ) = : b z" + b i z" ~ 1 + . . . 4 

i) Formula in terms of the traces: 

b =l ; b k = - l -{b k _ l \x(A)+b n ^r{A 2 ) + ...+b^(A k - x 

ii) Formula in terms of subdeterminants: 



-1Z- 



■b tx(A 



")} 



1 ; b k = (-lf £ det 

h<-<h 



a hh 



In particular, 



*JU2 



l }\]k 
2 hjk 



■ik I 



'1 



-tr(A) , b n = (-!)" det (A) 



V*=l, 



5 The signal cames from the number 1 + [k/2\ of lines of the matrix defining A* which are different by a factor —1 from the 
corresponding lines of the matrix defining Aj: by a basic property of the determinant, each different line implies a factor —1 in the 
relation between the two determinants. 
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Proofs for such formulas can be found in [7] and [8 1, respectively. Another formula can be found in (9)- 
Finally, I collect the facts in a simple algorithm: 

Algorithm 14 Characterization of positive operators 

Let A be a operator in C" and consider its representation in the canonical basis ( or any other orthogonal 
basis) of<C": 

W = ( a ij)i,j=l,...,n 

i) Verify if the operator is self-adjoint, i.e., ifdjj — aj, \/i,j — 1 , . . . , n; 

ii) Calculate the coefficients of the characteristic polynomial of A, using one of the formulas in Theorem \13\ 
Hi) Calculate the multiplicity of the number zero as a root of the characteristic polynomial of A, using, for 
example, the formula 

d k p A 



no = min < k G N / 



dz k 



^0 

=o 



iv) Calculate the Hurwitz determinants of the characteristic polynomial of A up to the order n — no, using 
Definition^ 

v) Verify if theses Hurwitz determinants satisfy the condition ofTheorem \12\ 

The computational efficience of this Algorithm [14] could be compared with the computational cost to 
calculate the operator's eigenvalues with sufficient precision. I let this for you. 



5 Special Cases 
5.1 Dimension n = 2 

Let A be a self-adjoint operator in C 2 and consider its matrix representation with respect to the canonical 
basis (or with respect to any other orthogonal basis): 

Self-adjointness means 

aij=dji e C , Vi,j = 1,2 

The coefficients of the characteristic polynomial of A are, as one can obtain directly from (fl]i or from Theo- 
rem[l3l 

b = 1 , b\ = -tr(A) , b 2 = det(A) 

The Hurwitz determinants (O are: 

Ai =b\ = -tr(A) , A 2 ^ bib i -b b 2 =det(A) 

Now, to get the necessary and sufficient conditions for A be positive, we have to consider three different situ- 
ations, distinguished by the multiplicity, denoted here by ji, of the number zero as a root of the characteristic 
polynomial of A. 

Case ji = (£>2 = det(A) ^ 0): TheoremfT2lstablishes that A is positive if, and only if, 

tr(A) >0 , det(A) >0 
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Case ji = 1 (b 2 = det(A) — , b\ — — tr(A) 7^ 0): Theorem[T2l stablishes that A is positive if, and only 

if, 

tr(A) >0 

Case ji = 2 (b 2 — det (A) — Q,bi — — tr(A) = 0): Theorem[T2lstablishes that A is positive (without any 
further condition). 

Finally, we can collect the cases in a single sentence: 

A self-adjoint operator A in C 2 is positive if, and only if, it satisfies one out of the two following 
conditions: 

f (0 det(A)=0,tr(A)>0 

\ (ii) det(A)>0, tr(A)>0 W 

5.2 Dimension n = 3 

Let A be a self-adjoint operator in C 3 and consider its matrix representation with respect to the canonical 
basis (or with respect to any other orthogonal basis): 

W = ("</ )»,;= 1,2,3 

As above, self-adjointness means 

a\j — dji^C , Vz, j = 1, 2,3 

The coefficients of the characteristic polynomial of A are, as one can obtain directly from (fl]i or from Theo- 
rem [T3] 

b = 1 

h = -tr(A) 

b 2 = -i(Z7 1 tr(A)+fe tr(A 2 )) = i(tr(A) 2 -tr(A 2 )) 

h = -^{Mr(A)+/7 1 tr(A 2 )+fe tr(A 3 )} = -det(A) 
The Hurwitz determinants of the characteristic polynomial are: 

A!=&i = -tr(A) 

A 2 =b 1 b 2 -b b 3 = ~2 tr ( A ) ( tr ( A ) 2 - tr ( A2 )) +det(A) 

A 3 = bib 2 h - b bl = j |tr (A) (tr (A) 2 - tr (A 2 ) ) - det (A) j det (A) 

Now, to get the necessary and sufficient conditions for A be positive, we have to consider four different 
situations, distinguished by the multiplicity of the number zero as a root of the characteristic polynomial of 
A. To save space (and also because it would be a little tedious), I consider only the case in which the number 
zero is not a root of the characteristic polynomial (bi = — det (A) ^ 0). In this case, Theorem[T2lstablishes 
that A is positive, if and only if, 

tr(A)>0 , itr(A) (tr(A) 2 +tr(A 2 )) +det(A) >0 , j^tr(A) (tr(A) 2 ~tr(A 2 )) - det (A) j det (A) > 
Acknowledgement 15 / thank my friend Wescley Bonomo for pointing out the Routh-Hurwitz criterion. 
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